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Abstract 

The quark viscosity in the quark gluon plasma is evaluated in HTL approximation. 
The different contributions to the viscosity arising from the various components 
of the quark spectral function are discussed. The calculation is extended to finite 
values of the chemical potential. 
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In this letter we present a microscopic calculation of the shear viscosity rj 
and of its ratio to the entropy density s for a system of quarks. This subject 
has been widely discussed in recent times, on the basis of the intriguing re- 
sults from the experiments carried out at the Relativistic Heavy Ion Collider 
(RHIC). In particular the measurement of the V2 coefficient [Tj in the mul- 
tipole analysis of the angular distribution of the produced hadrons seems to 
imply a very small viscosity, like the one of an almost perfect fluid [2|3lfi] . in 
contrast with the current description of the QGP as a gas of weakly interacting 
quasi-particles [5]. The substantial collective flow observed in these collisions 
also seems to imply quite small values for the viscosity prff8f9|T0] . 
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The temperature behavior of the ratio rj/s should determine the value Tc 
at which a phase transition occurs in the system [UTO]. Much emphasis has 
received the result of certain special supersymmetric Yang-Mills theories in 
4 dimensions, which predicts a lower limit for the viscosity/entropy density 
ratio, namely i]/s > 1/Att (in units h = ks = c = 1) fr\JT2\ : the latter value is 
approached by several recent calculations [T3|14|15p!6] . Transport coefficients 
from the linearized Boltzmann equation for a non abelian gauge theory have 
been obtained using perturbative QCD in Refs. [T7IIT8] . 

Here we shall describe the quarks in the QGP within the Hard Thermal Loop 
(HTL) approximation, which has already been proven quite successful for the 
thermodynamics of high temperature QGP [T^2Uf21j . Quark-antiquark corre- 
lators have also been widely considered in this approach (see, e.g., Refs. [22l[23p^ 
and references therein). 

We evaluate the shear viscosity t] of the system limiting ourselves to the quark 
degrees of freedom and starting from the Kubo formulas for the hydrodynamic 
transport coefficients; it is expressed in terms of the correlator of the energy- 
momentum tensor as follows: 
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where we consider the following retarded correlator (at zero momentum) [16j 
1 r d^h 
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In the above ivi = 27r///5, Un = (2n -|- l)7r//5, (/5 = 1/T) are Matsubara 
frequencies and the spectral representation of the massless quark propagator 
is employed 

S{tUJn,k) = - — , (3) 

J 2lT lUJn — OJ 
— oo 

the spectral function being conveniently split into the chirality projections: 

7° - y ■ k 7° + y • k 

P(^' ^) = n ^) ^ 7, P-(^' ^) • (4) 



The explicit expression of the HTL quark spectral function introduced in 
Eq.dl) reads p5]: 

p±(o;, k) = 27iZ±{k)[5{uj - lu±) + 5{uj + co^)] + 2np±{cu, k)e{k^ - lu^) , (5) 
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where nig = g(T)T/\^ is the quark thermal mass, 



Z±{k) 
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are the residues of the quasi-particle poles [see Eq. 
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The spectral functions p±{uj, k) consist of two pieces, a pole term and a cut: 

p^{u,k) = pT{u^,k) + pi'^iu.k) . (8) 

At a given value of the spatial momentum k, in the time-like domain (cu^ > k^) 
discrete poles are associated to quasiparticle (QP) excitations with dispersion 
relation u = uj±{k), while in the space-like domain {u"^ < k"^) a cut accounts 
for the Landau damping (LD) of a quark propagating in the thermal bath. The 
two poles correspond to quasi-particles with opposite chirality/helicity sign, 
their dispersion relation being given by the solutions of the implicit equation: 
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Vertex corrections to the correlator ([2]) deserve further investigation and will 
not be included here. 

We can now insert the spectral representation, Eq.Q, of the propagators 
into Eq.([2]) and sum over the Matsubara frequencies with a standard contour 
integration; then, by performing the usual analytic continuation iioi uj + it]^ 
(retarded boundary conditions), taking the imaginary part of the result and 
putting into Eq.([T]) we obtain [TSflG] : 
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where f{uj) = l/(exp'^'^'^ +1) is the thermal distribution for fermions with 
chemical potential p. 



The trace over spin can be explicitly carried out after inserting Eq. (jlj) into 
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Eq. (fTOj) and, since 



Tr 
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one finally gets: 



eOvr^T 

+ 00 



+00 



(12) 



J duj[l- f{uj)]f{uj) {pI{u, k) + p'_{uj, k) + 8p+{u, k)p4u, k)} 



By inserting the explicit expressions for p±{uj,k) into Eq. ( |T2l) this quantity 
turns out to be split, in the HTL approximation, into three different con- 
tributions [pole-pole term, labelled QP, pole-cut term (QPLD) which mixes 
quasi-particle and Landau damping parts and cut-cut (LD) which comes from 
Landau damping only]: 



The quasi-particle (pole-pole) contribution reads: 



-,QP . 
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dkk' J du[l-f{u)]f{u){lOp2''{oo,k)p^^{u,k)} .(14) 



The Landau damping (cut-cut) contribution reads: 



X 



607r3T J "^^^ J ^^[1-/H]/Mx (15) 

-k 

{[p^^(a;, k) - p^_^ico, k)f + lOpl^ico, k)p^^(c., k)} 
and finally the pole-cut contribution reads: 



4 



r/^^"^ = ^ J dkk'J [1 - /(.;)]/(..) X (16) 

-k 

X {2p?^(^, k)p^^(^, k) + 2p«^(a;, k)p^^(a7, k) 
+ 8p«^(a;, k)p^^(c., k) + 8p«^(c., k)p^^(a;, k)} . 

Among the three contributions, the QP one, as it stands, turns out to be 
divergent, as it corresponds to a free gas of particles. In order to cure this 
unphysical outcome, we shall assign, following Ref. [26], a finite width to 
the poles keeping unchanged the residues: the width of the poles ±u!+{k) is 
7+ and the width of ±uj^{k) is 7_. The spectral function for these dressed 
quasiparticles then reads: 

pfiuj,k) = Z±{k)-. - + Z^{kh ^J^i- ^ (17) 



and is normalized according to 

(18) 

The width is fixed as in Ref. 



:coramg to 

+ 00 

/ ^p2''iu;,k) = Z^{k) + Z_{k) 
s in Ref. ESI: 



7±=a%^. (19) 
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where a = 3.0 is chosen as an intermediate value between the one loop 
calculation and the improved one of Ref. [26] for Nj = 2. Cf = 4/3 is the 
Casimir constant for the fundamental representation of SU(iVj)Il] A word of 
caution is in order with the use of a constant quasi-particle width: indeed, as 
observed in Ref. [I6], the second term of rj^^ [second line in Eq. ( fT^ ] would 
contain an ultraviolet divergence. A natural cutoff is provided by the scale 
of momenta within the HTL approach, since the particles dressed by HTL 
self-energy are expected to be "soft" constituents, where the soft momenta 
are defined as the ones smaller than T (or gT). Hence the upper limit of the 
integration over k in the second line of Eq. (1141) has been set to T. 

In addition to the shear viscosity, we have then evaluated the entropy density 
according to the formulas of Ref. [20f2T] . 

The ratio rj/s is illustrated in Fig. [1] (solid line) for p = as a function of 
T/Tc, with Tc = 202 MeV, which was suggested as the transition temperature 



^ An additional reason for the present choice lies in the fact that in Ref. [26| the 
width is calculated at zero momentum, while we are integrating on the quasi-particle 
contribution over all momenta. 
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for Nf = 2 in Ref. [27]. Moreover we use the temperature dependent running 
gauge coupling given by the two- loop perturbative beta-function [5]|28ll29f23j . 
In the same figure we display the separate contributions from rj'^^, vp^^^ 
and rj^^ . Clearly the dominant contribution is 77*^^: we notice that the cutoff 
dependence in the second term of the QP contribution is quite mild, since 
the first term is by far dominant; indeed by artificially increasing the cutoff 
up to lOT the ratio rj/s increases by about 2%. We wish to stress that the 
contribution of the continuum to the viscosity was never evaluated in previous 
works: we find that it is quite small and that its temperature dependence is 
opposite to the QP one, showing a fast decrease, which is also reflected in the 
QP-LD interference term. 



00 




Fig. 1. The ratio r//s for = as a function of T jT^- the QP (dash-double-dotted 
line), LD (dashed) and the interference between QP and LD (dot-dashed) are shown, 
together with the total result (solid). The limit l/47r is indicated by the dotted line. 

We have then extended our calculation of s, 77 and 77/5 to the case of finite 
chemical potential, taking also into account the corresponding correction to 
the quark mass [30], = g^{T) (T^ -|- /i^/vr^) /6. Finite baryonic densities 

in high energy heavy ion collisions will be of interest for the forthcoming GSI 
experiments. The ratio of the entropy density to the corresponding Stefan- 
Boltzmann limit (sfj-ee) is shown in Fig. [2] for three different values of /i: the 
curves corresponding to /i 7^ have been extended to lower temperatures, 
since Tc should gradually decrease with increasing chemical potential. The 
effect of the chemical potential appears to be sizable only for the largest value 
considered here. 



In Fig. [3] we display the t]/s ratio for different values of /i: one can see that the 
effect of a finite chemical potential increases the ratio, which implies a larger 
correction to the viscosity than to the entropy density (the modification of the 
latter, according to the results shown in Fig. [21 should lower rj/s). Obviously 
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the largest effects are seen at lower temperatures. 
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Fig. 2. Entropy density of quarks in HTL approximation (divided by Sfree) as a 
function of T/Tc for /x = (solid line), ^ = 200 MeV (dashed) and /x = 400 MeV 
(dotted). 
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Fig. 3. The ratio ry/s as a function of T/Tc for different values of the chemical 
potential: /x = (sohd hue), [x = 200 MeV (dotted), n = 400 MeV (dot-dashed) 
and ^ = 600 MeV (dashed). 

In summary we have microscopically calculated the viscosity and the entropy 
density for a system of quarks described in the framework of the HTL approx- 
imation. The quasi-particle term in the spectral density has been modified by 
accounting, within the same approach, for a finite width, which removes the 
divergence in the shear viscosity derived from the Kubo formula. The results 
for the T]/ s ratio are comparable to the ones obtained within phenomenological 
models fT6] and turn out to be about a factor 3 larger than recent lattice cal- 
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culations [H] in the pure gauge sector. We found that the continuum part of 
the quark spectral density gives a neghgible contribution to the viscosity and 
the corresponding rj/s ratio slowly decreases with temperature, at variance 
with the dominant, quasi-particle one. 

We have also considered the effect of a finite chemical potential and we found 
that it causes a sizable increase, both in the entropy density and in the vis- 
cosity: indeed induces the largest corrections to the latter, so that also the 
ratio Tj/s tends to increase with the chemical potential. 
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